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Abstract 

Gravitational fields invariant for a 2-dimensional Lie algebra of 
Killing fields [X, Y] =Y, with Y of light type, are analyzed. The con- 
ditions for them to represent gravitational waves are verified and the 
definition of energy and polarization is addressed; realistic generating 
sources are described. 
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Introduction 



Gravitational waves, that is a propagating warpage of space time generated 
from compact concentrations of energy, like neutron stars and black holes, 
have not yet been detected directly, although their indirect influence has been 
seen and measured with great accuracy. Presently there are, worldwide, many 
efforts to detect gravitational radiation, not only because a direct confirma- 
tion of their existence is interesting per se but also because new insights on 
the nature of gravity and of the Universe itself could be gained. For these rea- 
sons exact solutions of the Einstein field equations deserve special attention 
when they are of propagative nature. The need of taking into full account 
the nonlinearity of Einstein's equations when studying the generation of grav- 
itational waves from strong sources is generally recognized 26J. Moreover, 
despite the great distance of the sources from Earth (where almost all the ex- 
perimental devices, laser interferometers and resonant antennas, are located) 
there are situations where the non linear effects cannot be neglected. This is 
the case when the source is a binary coalescence: indeed it has been shown 
jH] that a secondary wave, called the Christodoulou memory is generated via 
the non linearity of Einstein's field equations. The memory seems to be too 
weak to be detected from the present generation of interferometers j2Hl (even 
if its frequency is in the optimal band for the LIGO/ VIRGO interferometers) 
but of the same order as the linear effects related to the same source, thus 
stressing the relevance of the nonlinearity of the Einstein's equations also 
(soon) from an experimental point of view. 

On the theoretical side, starting from the seventy's new powerful mathe- 
matical methods have been developed to deal with nonlinear evolution equa- 
tions. For instance, a suitable generalization of the Inverse Scattering Trans- 
form allows to integrate jl] Einstein field equations for a metric of the form 

g — f (z, t) [alt 2 — dz 2 ) + h\\ (z, t) dx 2 + h 22 (z, t) dy 2 + 2hi 2 (z, t) dxdy. 

Indeed, the corresponding vacuum Einstein field equations reduce essentially 1 
to 

(aH^H^+fa H" 1 !!^ = 0, 

where H = \\h ab \\ , £ = (t + z)/y/2, rj = (t-z)/V2, a = </|detH|. 
This is a system of non-linear differential equations whose form is typical 

lr rhe function / can be obtained by quadratures in terms of the matrix H. 



1 



for two-dimensional integrable systems. Its solution through the Inverse 
Scattering Transform, yields gravitational solitary waves solutions. 

A geometric inspection of the metric above shows that it is invariant 
under translations along the x,y-a.xes, i.e. it admits two Killing fields, d x 
and d y , closing on an Abelian 2 two-dimensional Lie algebra Ai- Moreover, the 
distribution T>, generated by d x and d y , is 2-dimensional and the distribution 
T> orthogonal to T> is integrable and transversal to T>. 

Thus, it has been natural to consider the general problem of charac- 
terizing all gravitational fields g admitting a Lie algebra Q of Killing fields 
such that: 

I. the distribution V, generated by the vector fields of Q, is two-dimensional. 

II. the distribution T* 1 - orthogonal to T> is integrable and transversal to 

V. 

The aim of this article is to study, among those solutions, the ones which 
show a propagative nature. A preliminary account of our results is given in 



The article is organized as follows. In section 1 gravitational fields in- 
variant for a two dimensional Lie algebra are characterized by reducing the 
Einstein equations to the so called /i-deformed Laplace equation. Harmonic 
coordinates are also introduced. Section 2 is devoted to the analysis of the 
wave-like character of the solutions through the Zel'manov and the Pirani 
criteria. In section 3 realistic sources are described. In section 4 Landau- 
Lifshitz's and Bel's energy-momentum pseudo-tensors are introduced and a 
comparison with the linearised theory is performed. Eventually, section 5 is 
devoted to the analysis of the polarization of the waves. 

In the following, Kil (g) will denote the Lie algebra of all Killing fields of 
a metric g while Killing algebra will denote a sub-algebra of JCil (g). 

Moreover, an integral (two-dimensional) submanifold of T> will be called a 
Killing leaf, and an integral (two-dimensional) submanifold of V 1 - orthogonal 
leaf. 

2 The study of metrics invariant for a Abelian 2-dimensional Killing Lie algebra goes 
back to Einstein and Rosen jl()l I22| . Kompaneyets [121; recent analysis can be found in 

[23113 
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1 Invariant vacuum gravitational fields 



Let g be a metric on the space-time M and Q2 one of its Killing algebras 
whose generators X, Y satisfy the commutation relation 

[X,Y] = aY, a = 0,1. (1) 

The Frobenius distribution T> is two-dimensional and a coordinate system 
(x^ 1 ), /x = 1, 2, 3, 4 exists such that 

Such a coordinate system is called semiadapted (to the Killing fields). 

Condition II of the previous section allows to construct special semi- 
adapted charts such that the fields e\ = d/dx 1 , ei = d/dx 2 , belong to T> . 
In such a chart, called from now on adapted, the most general Q2- invariant 
metric has the form [23] 

g = g i jdx t dx^ + ^er 2 A (x 4 ) 2 — 2a fix 4 + isj dx 3 dx 3 + 
2 (/i — crAx 4 ) dx 3 dx 4 + Xdx 4 dx 4 , i,j = 1, 2 

where gij, A, /x, v are arbitrary functions of (x 1 ,^ 2 ). 

If the Killing field Y is not of light type, i.e. g(Y, Y) 7^ , condition II 
follows automatically from I. The local structure of this class of Einstein met- 
rics has been explicitly described; it turns out that a third Killing field 
Z exists such that (X, Y, Z) span the Killing algebra so (2, 1) and all invari- 
ant metrics are static and locally diffeomorphic to the pseudo-Schwarzschild 
metric. This metric was also found in the context of warped solutions [TT) l2"5]. 

In the following the Killing field Y will be assumed to be of light type. In 
this case the general solution of vacuum Einstein equations, in the adapted 
coordinates (x 1 , x 2 , x 3 , x A ), is given by 

g = 2f{dx 1 dx 1 + dx 2 dx 2 ) + fi[(w (x\ x 2 ) - 2ax 4 )dx 3 dx 3 + 2dx 3 dx% (2) 

where fi = -0$ + B; D, B G TZ, $ is a non constant harmonic function of 
x 1 and x 2 , f = (V$) 2 and w (x 1 ,^ 2 ) is a solution of the [i-deformed 

Laplace equation: 

Aw + (d x i In d x iw + (d x 2 In \fj,\) d x 2iu = 0. (3) 
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The solutions of the /z-deformed Laplace equation will also be called fi- 
harmonic functions. 

In the particular case a — 1, f — 1/2 and \l — 1, the above metrics 
are locally diffeomorphic to a subclass of the vacuum Peres solutions 
corresponding to a special choice of the harmonic function parameterising 
that metrics. 

This is easily seen by introducing, for u ^ 0, new coordinates (x,y,u,v) 
defined by x 1 = x, x 2 = y, x 3 = In \u\ , x A = uv. 

Of course, metrics (J2J) can also be obtained by imposing an additional 
isometry to metrics admitting a Killing vector field of light type and, indeed, 
they are an extension of the ones reported in ^ . 



When w is constant the above family of Einstein metrics admits the time- 
like Killing vector field 

T = exp (-x 3 ) A + I [( w - 2ax i ) exp (-a; 3 ) + exp (x 3 )] 

which is hypesurfaces orthogonal. This means that these Einstein metrics are 
just static gravitational fields. 

As it will be argued in the next sections, when w is not constant the above 
family of Einstein metrics may represent propagative gravitational fields. 

Since the distribution D 1 - is assumed to be transversal to £>, the restric- 
tion of g to any Killing leaf, say S, is non-degenerate. So, (S, g\ s ) is a ho- 
mogeneous two-dimensional Riemannian manifold. In particular, the Gauss 
curvature K (S) of the Killing leaves is constant. An explicit computation 
shows that K (S) vanishes. 

Thus, the space-time M. has a fiber bundle structure 

7i :M — >W, 

whose basis W is diffeomorphic to the orthogonal leaves and whose fibers are 
the Killing leaves and as such are flat two-dimensional Riemann manifolds. 

Despite the non-linear nature of general relativity, gravitational fields (0) 
obey to two superposition laws. Indeed, with two harmonic functions $i and 
$2 we can associate three gravitational fields (in facts a whole two-parameters 
family), that is, g<$ >1) g$ 2 and g a $ 1+ f,$ 2 ; the last one, which is associated with 
the linear combination of $i and $2, may be regarded as the superposition 
of the two associated solutions g§ 1 and g$ 2 . 

The second superposition law follows from the linearity of the /i-deformed 
Laplace equation, so that with two /i-harmonic functions W\ and W2 we 
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can associate three gravitational fields g Wl , g W2 and their sum g Wl+W2 = 
{gbun +92w 2 ) /2 

1.1 The harmonic coordinates. 

The coordinates (x 3 , x 4 ) on the Killing leaves S have a clear geometric mean- 
ing but are of difficult physical interpretation. Fortunately, being the Killing 
leaves flat manifolds, it is possible to introduce coordinates (z, t) diagonaliz- 
ing the metric 

g\s = V\{™ - 2ax A )dx 3 dx 3 + 2dx 3 dx 4 } , 

where fx and w, being the restriction of the functions fx and w to the Killing 
leaves, are constant. 

The coordinate system (x,y,z,t), where for z > t 

x = W\ (x , x 2 ) 
y = W2 (x 1 , x 2 ) 

z = i [{2x A — w (x 1 , x 2 )) exp (x 3 ) + exp (— x 3 )} 
t — 2 [(2x 4 — w (x 1 , x 2 )) exp (x 3 ) — exp (— x 3 )] , 

is harmonic, Wii^x 1 ^ 2 ) and w%{x 1 ,x 2 ) denoting any two independent /x- 
harmonic functions. The generic Killing leaf S is mapped onto the half-plane 
z > t, the line z — t representing the points with x 3 = +oo. 
In these coordinates, metrics (j2J) take the form 

g = 2 ^ i V$) J 2 UVy) 2 dx 2 + (Vx) 2 dy 2 - 2Vx-Vydxdy] (4) 
ix 

+H[dz 2 - dt 2 + d (w) d (In \z - t\})}, 

where J = d x iW\d x 2W2 — d x 2W\d x iW2 is the Jacobian determinant of the map 
(x 4 ,x 2 ) -> 

In the case fx = const, the /i-deformed Laplace equation reduces to the 
Laplace equation and W\, w 2 reduce to be just harmonic functions. Thus, 
it is possible to choose x = x 1 , y = x 2 so that in the harmonic coordinates 
(x,y,z,t), and for fx — 1, the above Einstein metrics take the particularly 
simple form 

g = 2f(dx 2 + dy 2 ) + dz 2 - dt 2 + d (w) d (In \z - t\\) . (5) 



5 



This coordinates system explicitly shows that, when w is constant, the Ein- 
stein metrics given by Eq. (j5J) are static and, under the further assumption 
$ = xy2, they reduce to the Minkowski one. Moreover, when w is not 
constant, gravitational fields (jSJ) look like a disturbance moving, along the z 
direction on the Killing leaves, at light velocity. However, the last observa- 
tion is neither rigorous nor covariant. Since the propagation direction is the 
most important ingredient in the study of the polarization, in the following 
sections a detailed analysis will be devoted to this question. 
In the following we will assume that w is not constant. 

2 Zelmanov's and Pirani's criteria 

To check the wave character of gravitational fields ©, a Zakharov gener- 
alization of the Zel'manov criterion 31 will be applied which states that 
a vacuum solution of the Einstein equations is a gravitational wave if the 
components R^xa of the corresponding Riemann tensor field R, satisfy a 
hyperbolic equation of the form 

g^VaVpR^xv = N^ Xa (6) 

where V/? denotes the Levi-Civita covariant derivative of the metric and 
NftuXa denote the components of a tensor field N depending at most on first 
derivatives of the Riemann tensor itself. 

For symmetric manifolds Eq. © with N = is an identity because the 
Riemann tensor is covariantly constant, but it may become an identity also in 
the case of Einstein manifolds {R a p = Kg a p), for special choices of the tensor 
field N. Hence, to exclude a priori these situations, the original Zel'manov 
criterion is formulated in the more restrictive assumptions: 

• RafijS not covariantly constant 3 ; 

• g a ^ a V p R pv xa = . 

The metrics in Eq. (J2J) certainly do not define symmetric or Einstein 
manifolds, as can be checked from the components of the Ricci tensor given 
below. Hence, the first hypothesis is certainly satisfied while the second 

3 That is the manifold is not symmetric. 
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one, i.e. N = 0, which ensures the applicability of the criterion to Einstein 
manifolds too, is not needed. 

Concerning the physical meaning of this criterion, it can be shown [31] 
that the characteristic hypersurface of the system of equations (0) is identical 
with the characteristic hypersurface of the Einstein and Maxwell equations 
in curved space-time. Consequently, Eqs. (JBJ) describe the propagation of the 
discontinuities of the second derivatives of the Riemann tensor. This links 
the Zel'manov criterion to the intuitive concept of local wave of curvature. 
The criterion is independent on the explicit form of N^xa] m fact, the char- 
acteristic hypersurface of a system of equations is determined only by the 
highest derivative term. Then we will not fix an explicit form of N^ u \ a but 
just require that N^xa be a tensor containing at most first derivatives of the 
Riemann tensor. This clearly corresponds to a covariant criterion [31] ■ Then 
a sufficient condition is 

g^dvdpR^ = 0, (7) 

where dp are the usual partial derivatives. In fact, if this is the case then 
N^vXa is a tensor containing at most first derivatives of the Riemann tensor. 
To start with, let us verify that gravitational fields (j3J) do satisfy Eqs. 

©• 

In the harmonic coordinates system the only nonvanishing components 
of the Riemann and Ricci tensor fields corresponding to metrics (0) are pro- 
portional to one of the following 



R 



txzx 



(2,fW, XX -\~ f \y W ,y flX^lX^) 



^f(z-t) 2 



(^ftVixy f ly l^ix f ix ) 

txzy ~ Af(z-t) 2 

(2fW m -f,yW,y+f jX W jX ) 

R 



tVZV ~ ^ 4/(*-f) a 

p fvy ^f^x ~f \fixx ~\~f?yy ) , g N 

r^xyxy ~J. j \Oj 



and 



respectively. 



Aw n (V/) 2 - fAf 



Rft 2 1 Rxx 



2f{z-tf 1 ^ 2f 
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Moreover, the harmonicity condition for $ implies the last component 
vanishes. In fact, when /x = const, A$ = implies for / that 

fAf - (V/) 2 = 0. (9) 

• When / is a constant function, Eqs. (JBJ) reduce to 

p _ W, XX _ W, X y _ W,yy , . 

tUxzx- 2(iZ _ f) 2> ^~2(z-t)2' 1 J 

which, w(x,y) being a harmonic function, are all harmonic functions 
of x, y. As a consequence, it is straightforward to check that the gen- 
eralized Zel'manov criterion, in the form (JJJ) , is satisfied [7j. 

• When / is not a constant function, the generalized Zel'manov criterion 
is still satisfied in the form thanks to a non trivial combination of 
the harmonicity condition for w and Eq. @. 

• In the general case when / and \x are not constant functions the Zel'manov 
criterion is satisfied in the form expressed by Eq. © • This may be more 
conveniently checked in the adapted coordinates x M . 

The Zel'manov criterion, even if it is covariant and allows a clear physical 
interpretation in terms of local waves of curvature, does not determine the 
propagation direction of the waves, that is the most important ingredient in 
the study of their polarization. In the next sections we will overcome this 
drawback of the Zel'manov criterion by using a suitable energy-momentum 
pseudo-tensor. 

Besides the Zel'manov-Zakharov criterion, the Pirani algebraic criterion, 
which is based on the Petrov classification, is satisfied. First of all, let us 
recall that a vacuum solution of the Einstein equations is a gravitational wave 
according to Pirani if its Riemann tensor is of type II, N or III in the Petrov 
classification JH] • Then, in light-cone coordinates (u = (z — t) / y/2, v = 
(z + t) I v2), where the metrics given by Eq.© read 

g = 2f(dx 2 + dy 2 ) + 2dudv + dw din \u\ , (11) 

the vector fields d u and d v are both isotropic. Moreover, it is trivial to show 
that the only non vanishing components of the Riemann tensor are 

Ruiuj i« o^ii^V 
ZU Z J 
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and this clearly corresponds to a type-N Riemann tensor in the Petrov classi- 
fication. Furthermore, it follows from the natural interpretation of the Pirani 
criterion [20] that the gravitational wave propagates along the null vector field 
d u , or, in other words, the gravitational wave (0) propagates along the z— axis 
with velocity c — 1. Thus, the Pirani criterion, even if with a less clear 
physical interpretation, allows an easy and covariant determination of the 
propagation direction. It will be an important self-consistency check for our 
calculations to discover the same results by means of the energy-momentum 
pseudo-tensors. 

3 The sources 

In the theory of gravitational waves a crucial problem is to characterise re- 
alistic sources able to generate waves enough strong to be detected by the 
experimental devices. 

The simplest source for metrics (J2J) (with a = 1) is dust with den- 
sity p and velocity U^ 1 and, then, characterized by an energy-momentum 
tensor T^ u = pU^U u . When is a light-like vector field, this kind of 
energy-momentum tensor can describe the energy and momentum of null 
electromagnetic waves, i.e., electromagnetic fields whose scalars, F^F^ and 
e a/3fiu F a pF^ v ^ both vanish. In this way it would be possible to describe the 
gravitational effects (in particular the emission of gravitational waves) of a 
very interesting astrophysical phenomenon, the 7— ray bursts (GRBs): emis- 
sion of ultra-high energetic 7-rays (~ 10 20 eV), whose origin is still to be 
fully understood. 

Moreover, with this source one can also preserve the symmetries [23 El 
of the vacuum solution. It can be shown that, in the adapted coordinates 
), the corresponding solutions of Einstein equations represent 
the following gravitational fields 

g = 2f(dx 1 dx 1 + dx 2 dx 2 ) + p[{w (x , x 2 ) — 2x 4 )dx 3 dx 3 + 2dx 3 dx 4 }, 

where p = D§ + B; D, B G 7Z, $ is a non constant harmonic function of x l 
and x 2 , f = (V$) 2 ^f\p\/ p, while w^x 1 ^ 2 ) is a solution of the p-deformed 
Poisson equation 

pAw + V/i ■ Vw = 2fp 2 p, 

in which, to save writing, beyond c = 1, it has been taken 8nG = 1, G being 
the Newton gravitational constant. 
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It is worth to remark that spin-1 gravitational waves in vacuum necessar- 
ily are not square integrable. Thus, if they are requested to be asymptotically 
flat, then 5-like singularities appear in the plane (x,y) in a neighborhood of 
the origin; in presence of matter sources the singularities are smoothed out. 

4 The energy-momentum pseudo-tensors 

The definition of momentum and energy associated with a gravitational field 
is an intrinsically controversial problem because these quantities are con- 
nected to the space-time translation invariance, whereas the group of invari- 
ance of general relativity is much bigger. With this cautionary remark in 
mind, various definitions are available which attain to different physical situ- 
ations. When dealing with the solutions of the linearised Einstein equations 
in the vacuum (plane gravitational waves) a commonly accepted definition is 
based on the canonical energy- momentum pseudo-tensor ( ED] ) : 



is the Ricci scalar deprived of terms containing the second derivative of the 
metric. 

Then, for the wave solutions of the linearised Einstein equations the en- 
ergy density Tq is expressed |0] as the sum of squares of derivatives of some 
metric components which do represent the physical degrees of freedom of the 
metric. 

Under a transformation preserving the propagation direction and the har- 
monic character of the coordinates system, in particular a rotation in the 
(x, y) plane, the physical components of the metric transform like a spin-2 
field. It is well known that in general in Eq. (|T2j) is not a tensor field 
but it does transform as a tensor field under those transformations which 
preserve the character of the field of consisting only of waves moving in the 
z direction, so that the g^ v remain functions of the single variable z — t. 

Thus, within the linearised theory, the canonical energy-momentum pseudo- 
tensor is a good tool to study the physical properties of the gravitational 
waves. 




dL 




where 




(13) 
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4.1 Comparison with the linearised theory 

The exact gravitational wave 

g = dx 2 + dy 2 + dz 2 -dt 2 + d (w) d (In \z - t\) , (14) 

given by Eq. (jSJ) for jj, — 1, / = 1/2 , has the physically interesting form of 
a perturbed Minkowski metric with h = dwdln \z — t\ . 

Moreover, besides being an exact solution of the Einstein equations, it is 
a solution of the linearised Einstein equations on a flat background too: 

f rTd li d v h = Q 

\ r)^{2h mv - h mp ) = 

Then, to study its energy and polarization, the standard tools of the 
linearised theory and in particular the canonical energy-momentum pseudo- 
tensor, could be used. Nevertheless, with h = d (w) d (In \ z — t\) the Tq com- 
ponent of the canonical energy-momentum tensor vanishes. This is due to 
the fact that the components of the tensor h cannot be expressed in the 
transverse-traceless gauge since h has only one index in the plane transversal 
to the propagation direction. 

Even if not explicitly declared, the standard textbooks analysis of the 
polarization is performed for the square integrable solutions of the wave- 
equation. Indeed, they can be always Fourier developed in terms of plane- 
wave functions with a light-like vector wave k^. 

The harmonicity condition for the plane wave solutions h^ v = e^e 1 ? + 
e *fiu e ~ ip with p = k^ and k^ = 0, reduces to 

~ = rTkue^x. (15) 

It is trivial to see that the symmetry group of this equation, which encodes 
the harmonic nature of the coordinate system, reduces to linear transforma- 
tions and more precisely to Poincare transformations [SHj (these are nothing 
but the usual "gauge transformations" of the linearised gravity). It can 
be easily shown that, for square integrable perturbations, one can always 
choose the transverse-traceless gauge. In other words, it is always possible 
to eliminate, with a suitable gauge transformation, the components of the 
perturbations with one index in the propagation direction (HOj, i-e. square in- 
tegrable perturbations of spin-1 do not exist. For these reasons, the canonical 
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energy-momentum pseudo-tensor, which is gauge invariant in the sense of the 
linearised gravity, "cannot see" the energy and momentum of gravitational 
fields given by Eq. (jl4|) which have one index in the propagation direction 4 
(but they are not square integrable so they cannot be gauged away). The 
fact that square integrable perturbations with one index in the propagation 
direction are always pure gauge is equivalent to the fact that, for such pertur- 
bations, the canonical energy-momentum pseudo-tensor identically vanishes. 



4.2 Landau— Lifshitz's and Bel's energy- momentum pseudo- 
tensors 

Besides the canonical energy-momentum pseudo-tensor, a deep physical sig- 
nificance can be given to the Landau-Lifshitz energy-momentum pseudo- 
tensor r pK ^3] defined by 

tPK = T^k < (2r ^ r - " rLr ^ " (/ V " " gPKgXn 

I pX pu f-pK per _j_ pK "per pre per pre per \ 

T i/ i/ \ L Act 1 fu> "T" 1 pu L Aer 1 ucr 1 Ap 1 Ap vol 

i kX pvlyp per | pp per _ pp per pP per \ 

^ i/ i/ V L Act 1 pi/ ^ 1 pi/ 1 Act 1 i/cr 1 Ap 1 Ap 1 va) 

+ g"Y v (rir; a -r p Xfl r K U(T )}. (ie) 

There are strong evidences that, in some cases, it gives the correct def- 
inition of energy [19] . In fact, the energy flux radiated at infinity for an 
asymptotically flat space-time, evaluated with the Landau-Lifshitz energy- 
momentum pseudo-tensor, has been seen to agree with the Bondi flux [H] 
that is with the energy flux evaluated in the exact theory. 

It is easy to check that the components p p = Tq of the 4-momentum 
density are 

P° = (74* K-) 2 + C ^v) 2 ] + (I^j^V-tlVH 2 V4 
p 1 = p 2 = 0, p 3 = p°, 

where Cj are some positive numerical constants, V = (d x ,d y ) and the har- 
monicity condition for w has been used. 



4 Of course, it is possible to find |23j a coordinate system in which the perturbation h 
has non vanishing components only in the transverse plane. However such a coordinate 
system is not harmonic. 
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The use of the Bel's superenergy tensor j3] 

rpaf3Xfi ^_ ^j^apXcr j^f3 p _|_ * j^apXa * p ^ 

where the symbol * denotes the volume dual, leads to the same result. Indeed, 
in adapted coordinates the metric has the form 

g = dx 1 dx 1 + dx 2 dx 2 + (w (x, x 2 ) — 2x 4 )dx 3 dx 3 + 2dx 3 dx 4 

and the only non vanishing independent components of the covariant Rie- 
mann tensor R a f3~/6 = 9a P R P p^s are 

-R1313 = —w,i\ ; -R1323 = —w,\2 ; -R2323 = —w,22 ■ 

It follows that the density energy represented by the Bel's scalar 

W = T al3Xfl U a U^U x U^ 

the U a 's denoting the components of a time-like unit vector field, depends on 
the squares of w,ij. Thus, both the Landau-Lifshitz pseudo tensor and the 
Bel superenergy tensor single out the same physical degrees of freedom. In 
particular, we can take the components ht x and ht y as fundamental degrees 
of freedom for the gravitational wave (|14j). 

Concerning the definition of the polarization, the above form for Tq is par- 
ticularly appealing because, apart from a physically irrelevant total derivative 
that does not contribute to the total energy flux, the component Tq repre- 
senting the energy density is expressed as the sum of square amplitudes. The 
momentum p l = Tq is non vanishing only in the z— direction and it is propor- 
tional to the energy with proportionality constant c = 1; that is these waves 
move with light velocity along the z— axis. Moreover, this result is perfectly 
consistent with the one obtained with the Pirani criterion. 



5 Spin 

The definition of spin or polarization for a theory, such as general relativity, 
which is non-linear and possesses a much bigger invariance than just the 
Poincare one, deserves a careful analysis [21] . 

It is well known that the concept of particle together with its degrees of 
freedom like the spin may be only introduced for linear theories (for example 
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for the Yang-Mills theories, which are non linear, it is necessary to perform a 
perturbative expansion around the linearised theory) . In these theories, when 
Poincare invariant, the particles are classified in terms of the eigenvalues 
of two Casimir operators of the Poincare group, P 2 and W 2 where P M are 
the translation generators and = \e^ vpa P v M pa is the Pauli-Ljubanski 
polarization vector with M^ v Lorentz generators. 

As it has been shown, the gravitational fields of Eq. (|14j) represent grav- 
itational waves moving at the velocity of light, that is, in the would be 
quantised theory, particles with zero rest mass. Thus, if a classification in 
terms of Poincare group invariants could be performed, these waves would 
belong to the class of unitary (infinite-dimensional) representations of the 
Poincare group characterized by P 2 = 0, W 2 = 0. But, in order for such 
a classification to be meaningful, P 2 and W 2 have to be invariants of the 
theory. This is not the case for general relativity, unless we restrict to a 
subset of transformations selected for example by some physical criterion or 
by experimental constraints. For the solutions of the linearised vacuum Ein- 
stein equations the choice of the harmonic gauge does the job [30J. There, 
the residual gauge freedom corresponds to the sole Lorentz transformations. 

For these reasons, only gravitational fields represented by Eq. (JTljl will be 
considered, which, besides being exact solutions, solve the linearised vacuum 
Einstein equations as well. There exist several equivalent procedures to eval- 
uate their polarization. For instance, one can look at the Tq component of 
the Landau-Lifshitz pseudo-tensor and see how the metric components that 
appear in Tq transform under an infinitesimal rotation 1Z in the plane (x, y) 
transverse 5 to the propagation direction 6 . 

The physical components of the metric are h tx and h ty and under the 
infinitesimal rotation TZ in the plane (x, y) transform as a vector. Applied to 
any vector (v\, Vq) the infinitesimal rotation 1Z, has the effect 

TZvi = v 2 , Uv 2 = -t>i, 

from which 

TZ 2 Vi = -Vi i = 1, 2, 
so that ilZ has the eigenvalues ±1. 

5 With respect to the Minkowskian background metric, this plane is orthogonal to the 
propagation direction. With respect to the full metric this plane is transversal to the 
propagation direction and orthogonal only in the limit \z — 1\ i— > oo. 

6 It has been said before, that this transformation preserves the harmonicity condition. 
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Thus, the components of h^ v that contribute to the energy correspond to 
spin-1 fields 7 , provided that only Lorentz transformations are allowed. 

Spin-0 and spin-1 gravitons have been considered, in a different context, 

in 0d CHI 123- 

6 Conclusions 

It has been shown that gravitational fields (jllj) represent spin-1 gravitational 
waves and that the reason why it is commonly believed that spin-1 gravi- 
tational waves do not exist is that, in dealing with the linearised Einstein 
theory, all authors implicitly assume a square integrable perturbation. In 
other words, square integrable spin-1 gravitational waves are always pure 
gauge. However, it has been proven that there exist interesting non square 
integrable wave-like solutions of linearised Einstein equations that have spin- 
1. These solutions are very interesting at least for two reasons. Firstly, they 
are asymptotically flat (with at least a 5-like singularity) in the plane or- 
thogonal to the propagation direction. Secondly, they are solutions of the 
exact equations too, so that the spin-1 cannot be considered as an "artifact" 
of the linearised theory. Realistic sources able to smooth out the mentioned 
singularities have also been found. 
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